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Abstract 

The dynamical properties of exciton transfer coupled to polariza- 
tion vibrations in a two site system are investigated in detail. A 
I fixed point analysis of the full system of Bloch - oscillator equa- 

tions representing the coupled excitonic - vibronic flow is per- 
formed. For overcritical polarization a bifurcation converting the 
stable bonding ground state to a hyperbolic unstable state which 
is basic to the dynamical properties of the model is obtained. The 
phase space of the system is generally of a mixed type: Above 
bifurcation chaos develops starting from the region of the hyper- 
bolic state and spreading with increasing energy over the Bloch 
sphere leaving only islands of regular dynamics. The behaviour 
pj of the polarization oscillator accordingly changes from regular to 

chaotic. 

O PACS number(s): 05.45. + b, 63.20.Ls, 82.20.Rp 
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1 Introduction 



The coupling between electronic and vibronic degrees of freedom in molecu- 
lar and condensed media is one of the basic mechanisms influencing transfer 
properties of electronic excitations in these systems. The investigation of its 
consequences started from the polaron problem (see e. g. [1] and references 
therein; for exciton - phonon interaction see [2]) and continued with the study 
of the influence of the vibronic bath variables on the excitation transfer prop- 
erties in the framework of the generalized Master equation [3] and stochastic 
Liouville equation approaches [4]. 

With the development of the theory of dynamical systems in the last decade 
(see e. g. [5], [6]) it has become attractive to analyze the implications of this 
basic mechanism employing concepts and methods of this latter field. Using a 
dynamic system approach one focuses on the detailed picture of the real time 
evolution of the relevant variables of the system. Of particular importance 
in this connection is the by now well established fact that systems generated 
by Hamiltonians with coupled degrees of freedom in many cases possess state 
spaces of a mixed type: regions of regular motion are embedded in a chaotic 
environment (see e. g. [7]). This results in completely different dynamic be- 
haviour of the system variables depending on initial conditions and / or system 
parameters. 

Employing such a point of view in this paper the dynamics of one of the 
simplest systems with coupled electronic and vibronic degrees of freedom is 
analyzed, namely the dissipationless transfer of an electronic excitation between 
two sites when each site is coupled to a vibration representing the polarization of 
the environment. To be definite we will consider the transfer of a Frenkel exciton 
between the two monomers of a molecular dimer where in each of the monomers 
the exciton is coupled to an intramolecular vibronic degree of freedom. The 
Frenkel exciton can of course be substituted by any other quasiparticle described 
by the same Hamiltonian. In extension of earlier approaches to the problem, 
where the effect of the vibronic degrees was modelled by a cubic nonlinearity 
(the DST-dimer originating with the paper [8] - the bifurcation of the flow 
line picture into a homoclinic structure and the appearance of chaos for the 
perturbed DST dimer were considered in [9], [10]) we will take full account of the 
dynamics of the vibronic variables, too. Hence our results demonstrate how the 
two state DST results are "embedded" in a more complete dynamic approach. 
For the effect of nonlinear self - trapping this problem of how DST results are 
embedded in a complete dynamical approach is of considerable interest. As an 
example we quote self - trapping at nonlinear impurity oscillators in a linear 
host [11]. 

The coupled excitonic - vibronic system is treated in a mixed quantum - 
classical approximation. As has become evident recently the coupling of clas- 
sical and quantum degrees of freedom can lead to chaos in both the quantum 
and the classical subsystems even in the simplest models of this type [12]. The 
mixed quantum - classical description is obviously justified whenever the quan- 
tum effects in one subsystem are negligible compared to the other. Such a 
description, however, can also constitute an important step towards the quanti- 
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zation of the whole system when the full quantization poses a severe challenge 
and the classical treatment of one part has to serve as a necessary guide. The 
most prominent example of an approach of this kind is the well known Born - 
Oppenheimer approximation. Having established a state space of both, regular 
and chaotic dynamics for the system in such a mixed description in a next step 
all the methodology of the "quantum signatures of chaos" [13] can be applied 
to the full quantization of the system. A system for which the correspondence 
between non - integrability in the mixed quantum - classical description and the 
properties of the fully quantized version has been already studied is e. g. the 
two - level atom coupled to an electromagnetic field with classical and quan- 
tum treatment of the latter (see e. g. [14]). The treatment of the excitonic 
- vibronic system in the quantum - classical approximation in our paper com- 
pletely confirms the conclusions drawn in [12] that the simplest systems in this 
approximation can display chaos in both subsystems and furthermore provides 
an example for the field of excitonic - vibronic coupled systems. 

In section 2 the model and the basic equations are formulated. Section 3 
contains the fixed point analysis and the related bifurcation. In section 4 we 
present an analysis of the dynamic properties using Poincare sections of the 
coupled excitonic - vibronic flow. 



2 Model and basic equations 

We consider the full dynamics of a coupled exciton - vibronic system where the 
excitaton is moving between the two monomers of a molecular dimer and at 
each monomer the exciton is coupled to an intramolecular vibronic degree of 
freedom. Our model is specified by the following Hamiltonian 

H = H exc + H v + Hi nt , (1) 

where H exc , H v and H v n are the excitonic, vibronic and interaction parts, 
respectively. H exc represents a standard two site model 

n n,m 

where c n is the probability amplitude of the exciton to occupy the n-th molecule 
and V nm the transfer matrix element due to standard dipole-dipole interaction, 
here and in what follows n, to = 1,2. The vibrational part H v is taken as the 
sum of the energies corresponding to intramolecular vibrations at each of the 
monomers for which we use the harmonic approximation 

H v = Y.\(A + ^lti)- (3) 

n 

Here q n ,p n and uj n are the coordinate, the canonic conjugate momentum and 
frequency of the intramolecular vibration of the n-th molecule, respectively. The 
interaction Hamiltonian takes into account that the exciton energy depends on 
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the molecular configuration of the monomers as expressed by the coordinates 
q n . Using the corresponding expansion up to first order in q n one gets 

Hint = ^2,-lnqnC* n C n , . (4) 

n 

where the j n are the coupling constants. In the case when several intramolecular 
modes are coupled to the exciton an additional summation has to be introduced 
in (3), (4). Being interested in the principal effects of the excitonic - vibronic 
interaction we restrict the consideration to the coupling with one configurational 
coordinate at each molecule. In what follows we will use a mixed quantum 
- classical description of the dynamics of the system, describing the exciton 
in the quantum context and the vibronic degrees of freedom in the classical 
approximation. 

From (1) - (4) one obtains the coupled equations of motion 



tC r , 



{^n ln(ln) c n ~\~ ^ ^ Vnm c m (5) 



4n = Pn (6) 
Pn = ~^ 2 q n ~ ln\c n \ 2 (7) 

(we set h = 1). 

One way to simplify the dynamics following from the system (5 - 7) is to 
assume quasistationarity in the vibrating subsystem which leads to the DST 
approximation mentioned in the introduction. Formally this is achieved by as- 
suming for the vibronic part q n = 0, p n = which results in q n = — (j n /uj 2 )\c n \ 2 
(the vibronic coordinate q n is "slaved" by the excitonic occupation |c n | 2 ). In- 
serting the latter expression into (5) one obtains 



ic, n 



\n\^n\ )^ri~t~ ^ ^ V~n 



where Xn will be given in (16). Eq. (8) is called the discrete self trapping (DST) 
equation. The transfer dynamics resulting from it were intensively analyzed for 
different parameter sets and numbers of molecular units N,n = 1,..N (see e. 
g. [16]). An equation of the type of (8) appears already in the early papers 
on the polaron problem as a limiting case [1]. One of the aims of this paper 
is to relate the dynamics following from eq. (8) for the dimer, i. e. N = 2, 
which was studied in [9], to the dynamics of the full system (5) - (7). We note 
that an intuitive statement concerning the validity of the DST approximation 
is usually given in the form that it corresponds to the case of fast polarization 
/ slow transfer and that then the vibronic subsystem is always in its ground 
state. One way to achieve this is the introduction of a fast relaxation into the 
vibronic subsystem. This, however, is beyond the scope of this paper where we 
remain on the dissipationless Hamiltonian level of description. For the role of 
dissipation in connection with the DST approximation we refer to the recent 
paper [15] and the references therein. 

Introducing the density matrix elements 

Pmn — ^n^ni (9) 
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and passing to the Bloch variables 



P12 + P21 

i(P2i~Pi2) (10) 
P22 ~ P11 



one obtains from (5) to (7) for the excitonic subsystem the Bloch equations 

x = — Ae(t)y 

y = z + Ae(t)x (11) 



z = -y 



where Ae(t) is given by 



Ae(t) := 72tt(*)-7igi(*) + *2-6i (12) 

In (11) we have introduced a dimensionless time by setting t 2Vt. Obviously, 
the dimensionless quantity Ae(t) controlling the coupling between the excitonic 
and vibronic subsystems is the dynamic generalization of the difference between 
the excitonic energies at each of the molecules caused by the interaction with 
the intramolecular oscillators. For a two site system (11) is equivalent to (5) 
but has the advantage of excluding the global phase. 

The coupled system of Bloch - oscillator equations (6), (7), (11) has two 
integrals of the motion: the energy H = E and the radius of the Bloch sphere 
x 2 + y 2 + z 2 = 1. Correspondingly, the number of independent variables of 
the system is reduced by two resulting in a total of five. A further reduction 
of the number of variables can be achieved by observing that only an effective 
oscillator, which is a linear combination of q\ and q^t enters the expression for 
Ae(t) : After the introduction of dimensionless oscillator coordinates Q n : = 
^q n , momenta P n := ^p n and a dimensionless energy difference 

a:=^il (13) 

the quantity Ae(t) is representable in the form 

Ae(t) = Q + a (14) 

where 

Q ■= X2Q2 - xiQi (15) 

and 

Xn ■= lll^l- (16) 

were introduced. By using the oscillator equations (6) and (7) it is easily verified 
that the equations for Q and the corresponding momentum 

P ■= X2P2 ~ X1P1 (17) 
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close, if the additional condition u\ = LO2 is satisfied. In what follows we 
will assume that this condition is fulfilled, i. e. the oscillators have identical 
frequencies. We stress that the dimer can still be "electronically disordered" 
when it has different excitonic energies e n and coupling constants j n . Using 
the condition u\ = loi one obtains the equations of motion for the effective 
oscillator 

Q = rP 

P = - r Q- r (pz + q) (18) 



where the dimensionless paramters 



P 



X2 + Xi 

w 



q := (19) 
r := lo/2V 

were introduced, r is a parameter translating between excitonic and vibronic 
time scales. Assuming quasistationarity for the effective oscillator and inserting 
the resulting equation Q = — (pz + q) into (14) one obtains the Bloch subsys- 
tem (II) in DST approximation as in [9]. The five equations (11) and (18) 
coupled via (14) still have two integrals of the motion: besides the radius of the 
Bloch sphere the dimensionless energy e := E /V of the considered subsystem 
is conserved. It can be expressed as the sum of the contribution from an effec- 
tive oscillator, whose equilibrium position depends on the current state of the 
excitonic subsystem, and a second exclusively excitonic term: 

e = - x - ( q - a )z - ^ + ^ \(Q + q + pz) 2 + P 2 ] (20) 

We note, that the condition imposed on Q and P to yield the DST approxima- 
tion implies that the effective oscillator has zero energy. For the DST dimer it 
is therefore straightforward to obtain from (20) the analytical solution for the 
orbits on the Bloch sphere. For the full dynamics which we consider here this 
can be generalized to an energy dependent restriction for the accessible part of 
the Bloch sphere: Since the energy of the effective oscillator is non-negative, 
the total energy is an upper bound to the excitonic term. This results in 

x > ~7^ z2 ~ (l ~ a ) z ~ e i (21) 
where the equal sign just describes the DST orbits. 



3 Fixed point analysis and bifurcation 

3.1 General equations 

We perform the fixed point analysis including the case of a nonsymmetric dimer 
(e\ ^ e 25 7i 7^ 72)- Setting in the equations of motion the time derivatives of z, 



6 



Q and P equal to zero, we find for the stationary states 

y s = o 

Q s = -(pz s + q) (22) 
P.s = 

The values for x s and z s are determined by y = 0. Using (22) one obtains 

(a - q)x s + z s - px s z s = (23) 

The restriction of the excitonic variables to the Bloch sphere yields the addi- 
tional condition 

x 2 s + z 2 s = l. (24) 

In order to determine the character of the fixed points we need the lineariza- 
tion of the dynamics for small deviations from them. The eigenvalues of the 
corresponding stability matrix turn out to obey the characteristic equation 

A 4 + (a?; 2 + r 2 )A 2 + r 2 (x; 2 - px s ) = 0. (25) 

with the location of the fixed point still general. According to the Hamiltonian 
character of the problem it is biquadratic and the sum of all stability exponents 
A is zero. They come in pairs 

A 2 = _ x ° 2 + r2 ±ly/(x7 2 -r*y + 4pr*x s . (26) 

It is appropriate to subdivide all stationary points according to whether they 
are located in the bonding (x s > 0, containing the bonding state x s = +1) or 
antibonding region (x s < 0, containing the antibonding state x s = —1). There 
is no transition between the two groups since x s = is excluded by (23). 

Let's consider the bonding region first. From (26) is evident that in this 
region A 2 is always real for both pairs of stability exponents and negative for 
at least one of them. The sign of the other one can be positive or negative. In 
the latter case all stability exponents are purely imaginary and the fixed point 
is stable elliptic. The condition for this to hold is 

P < ^~ 3 - (27) 

If it is not fulfilled the fixed point has a hyperbolic manifold and is unstable. 
In the antibonding region A 2 itself can be complex if 

\r 2 - x~ 2 \ < 2r^p\x s \. (28) 

In this case the antibonding fixed point will be unstable. 

From (22) it is clear that the effective oscillator in (20) has zero energy 
at every fixed point. Consequently, the fixed points of the full problem are 
identical with those from the DST problem [9] and we can use the DST orbits 
from (21) to show their location on the Bloch sphere. We stress, that the 
equivalence to the DST approximation does not hold for the stability: The 
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applicability of a time scale separation leading to the DST dimer depends on 
the value of the parameter r which does enter neither the location of fixed points 
(23) nor the stability condition for the bonding region (27) but the stability in 
the antibonding region (28). Consequently, the fixed point in the antibonding 
region was found to be always stable in [9] whereas it can be unstable in our 
more general approach. 

3.2 Symmetric dimer 

Finding the solutions of (23) it is convenient to consider the case of a symmetric 
dimer (e x = e 2 , Xi = X2 5 i- e - a = q = 0) first. Then the stationarity condition 
(23) depends on the parameter p only and reduces to the form 

(l-px s )z s = (29) 

We note that (29) also includes the situation when there is an exact compen- 
sation of the energy mismatch a and the polarization difference q. Considering 
the solutions of (29) one has to take account of (24) restricting x s and z s to 
values x s < 1 and z s < 1. So one finds the following solutions in dependence 
on the value of the parameter p: 

- Bonding region (x s > 0): 

• Case A: < p < 1 In this case the only solution of (29) is z s = and 
one finds from (24) the stationary point g with x s = +1. The point g is 
the ground state corresponding to a symmetric combination of excitonic 
amplitudes c\ = c 2 = l/\/2- The stability condition (27) is fulfilled. 

• Case B: p = 1 In this special situation we find two additional fixed points 
from the solution of (29) which are degenerate with the point g. 

• Case C: p > 1 There are two solutions of (29): z s = and x s = 1/p. 
From (24) we get for the first case x s = ±1 and for the second case 
z s = ±-\/l — (1/p) 2 =: ±^o- Correspondingly, we obtain three stationary 
points 

gl: x s = 1/p; z + = +z 
g2: x s = 1/p; z~ = -z 
h: x s = +1; z s = 

The points gl and g2 are stable, whereas h is unstable hyperbolic. 

The comparison of the structure of stationary points in the cases A and C 
reveals that the parameter p governs a pitchfork bifurcation: The ground state 
g below the bifurcation (p < 1) splits into two degenerate ground states gl and 
g2 above bifurcation (p > 1). At the former ground state a hyperbolic point h 
appears. This is also evident from fig. 1. The DST trajectories shown there can 
be interpreted in the full problem as isoenergy level lines within the subspace 
of zero oscillator energy. The physical reason for the bifurcation is that for 
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p > 1 the polarization is strong enough to hold the excitation preferentially at 
one of the molecules of the dimer. This corresponds in fig. 1 (b) to closed lines 
around the degenerate ground states. The regions on the Bloch sphere encircled 
by them can never be left by the exciton since the accessible part of the Bloch 
sphere as given by (21) will shrink further if the vibronic subsystem is not in 
the state of minimal energy. In this way we can conclude from fig. 1 (b) that 
above bifurcation there exist trajectories in the higher dimensional phase space 
of the full excitonic - vibronic problem where the exciton is self trapped, i. e. 
with a time average of z different from zero. 

- Antibonding region {x s < 0): In this region the only solution of (23) 
is z s = 0, x s = 1. This is the excited stationary state e and corresponds 
to an antisymmetric combination of excitonic amplitudes c\ = —c^ = l/y/2. 
According to (28) this point is stable for ^ r j 1 > 2^/p which will hold for 
the limiting cases r < 1 and r 1. If, however, the phonon and transfer 
frequencies are close to each other (r ~ 1) the stability is lost. 

3.3 Asymmetric dimer 

We now turn to the solution of eq. (23) for the general case of an asymmetric 
dimer. Without loss of generality the difference q — a can be considered one 
parameter. Therefore we will absorb the parameter a into a new defined q : = 
q — a. Inserting x s = ±^1 — z 2 s into (23) one obtains the following cases: 

- Bonding region {x s > 0): 

q + pz s = , Zs (30) 
V 1 " Z s 

In this case there are one and three solutions possible. In fig. 2 we present plots 
corresponding to these cases. In between there is the special situation when the 
straight line on the l.h.s. of (30) is tangent to the r.h.s. The geometrical 
condition of tangency in the intermediate situation can be used to derive the 
bifurcation diagramm fig. 3 in the q, p parameter plane. One obtains two lines 
±g(p) separating the region with a single solution for (30) from the one with 
three solutions. The lines are given by 

q(p) = (p-p 1/3 )(l-p- 2/3 ) 1/2 , P>1 (31) 

Starting from the symmetric dimer and using (27) one can show that in the 
case of three fixed points two are always stable and one is unstable hyperbolic. 
If there is just one fixed point, it must be stable elliptic. Thus the bifurcation 
diagramm generalizes the bifurcation of the ground state as expressed by the 
three cases A, B and C for the symmetric dimer to the asymmetric case. In 
particular, for p > 1 and q within the area restricted by the two lines in fig. 3 
the former ground state has split into two new elliptic and one hyperbolic state. 
In the asymmetric case these elliptic states have, however, different energies, 
i. e. the ground state degeneracy is destroyed. We note the interesting effect 
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that with increasing asymmetry, i. e. increasing value of \q\, the region of 
three fixed points is left and consequently the bifurcation is inverted. This is 
also evident from fig. 4 where the positions of z s and x s for a fixed p > 1 
are presented in dependence on q. We note that for a large linear asymmetry 
realized by \a\ \q\ ~ p the solutions of (23) approach the values for the linear 
asymmetric dimer (x s — ► ±(1 + a 2 ) -1 / 2 — ► 0, z s — ► ±a(l + a 2 ) -1 / 2 — ► ±1). 

- Antibonding region (x s < 0): 



This equation has only one solution corresponding to the excited stationary 
state. It is stable or unstable according to whether (28) is fulfilled or not. 

The location of the fixed points and the bifurcation for the asymmetric case 
is illustrated in fig. 5 for one special value of the asymmetry parameter q. 

As already stated above, the stationarity condition for the full set of vari- 
ables contains DST stationarity as a special case. Comparing the bifurcation 
with the one obtained in the DST approximation of the transfer problem (see 
[9]) one therefore finds that it is of the same pitchfork type represented by a 
ground state splitting and the appearance of a hyperbolic point. 

4 Analysis of the phase space flow 
4.1 General remarks 

Before turning to a detailed investigation of the dynamical regimes of our model 
we mention some general points concerning the coupled excitonic - vibronic flow. 
First of all we note that this flow is realized on a foliated phase space consisting 
of the surface of the Bloch sphere and the phase plane representing the dynamic 
states of the exciton and the oscillator, respectively. The flow of the excitonic 
and vibronic variables is restricted to each of the subspaces which constitute 
the "leaves" of the foliated space. On each of these two dimensional subspaces 
the dynamics of a one degree of freedom subsystem is realized. Taking account 
of the energy conservation three of the four dynamic variables of the coupled 
systems are independent. Hence by fixing one of the three independent variables 
a Poincare section on a two dimensional submanifold is determined. Depending 
on the choice of the independent variables and which of them is fixed we have 
Poincare sections realized on (i) the surface of the Bloch sphere, (ii) the phase 
plane of the oscillator and (iii) mixed manifolds. 

In the next subsection we demonstrate the characteristic dynamic regimes 
connected with the presence of the bifurcation and in particular the appearance 
of a hyperbolic point at the former ground state. These effects are displayed by 
the symmetric as well as by the asymmetric dimer and hence are presented for 
the simpler symmetric case. For the same reason we demonstrate the effect of 
the resonance in the simplest case of a 1-1 resonance for the symmetric dimer. 



q + pz s = - 
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Considering the asymmetric dimer we found flow pictures which may be 
thought as a continous deformation of the symmetric case. They do not con- 
tain qualitatively new results and are therefore not displayed. However, there 
are some special effects connected with the two dimensional character of the 
bifurcation diagram which will be discussed in the last subsection. 

4.2 Symmetric dimer 
Below bifurcation 

In fig. 6 (a) and (b) Poincare sections on the surface of the Bloch sphere are 
presented for two different energies and below the bifurcation (p = 0.5 < 1). 
They correspond to the left turning point of the effective oscillator (P = 0, 
(IP I dt > 0). Poincare sections at the right turning point contain the same 
information as one can see from a symmetry property: From every trajectory 
in the symmetric dimer another trajectory can be found by replacing z and Q 
with — z and —Q, respectively, since the total energy (20) is invariant under 
this transformation. From the Poincare sections one finds regular dynamics for 
both, low and high energy. For low energy the elliptic structure is centered in 
the region of the bonding ground state (x = 1, (f> = 0) whereas for high energy 
an analogous elliptic structure in the the region of the antibonding excited state 
(x = —l,(f> = 7r) is obtained. The oscillator also displays regular dynamics as 
is shown in the Poincare section of fig. 7. 

Above bifurcation 

In figs. 8 (a) - (d) a typical set of Poincare sections corresponding as in fig. 6 
to the left turning point of the oscillator but for a polarization strength above 
bifurcation (p = 2.0 > 1) is shown. For low energy one finds regular dynamics 
in the region of the new ground states (x s = l/p,z = ±^o) as m fig- 6a. 
These regular trajectories are "trapped" solutions of the excitonic subsystem 
in which the exciton is preferentially at one of the sites of the dimer. Their 
time average of z is 1 ~ ^fzo for the exciton being trapped at the first or second 
site, respectively. Correspondingly, the polarization oscillator performs small 
oscillations around the equilibrium positions Q = ±pz$. 

For intermediate energies involving the hyperbolic point, i. e. the dimen- 
sionless energy e = — 1 which can be used as a reference energy separating the 
low and high energy regions, one still obtains regular behaviour. From fig. 8 
(b) the structure of the hyperbolic point is evident. In fig. 8 (c) and (d) the 
spreading of chaotic dynamics over the Bloch sphere with increasing energy is 
shown. Beginning at the hyperblic point more and more trajectories become 
detrapped, i. e. the time average for z tends to zero with increasing energy. 
Chaos starts in the vicinity of the hyperbolic point (x = 1, z = 4> = 0, e = — 1) 
and then spreads over the Bloch sphere leaving only regular islands in the region 
of antibonding states. We observe that for chaotic trajectories the energy of the 
polarization oscillator is relatively high whereas for the regular orbits around 
the antibonding states the energy of the excitonic subsystem is at its highest 
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level and the oscillator does not get enough energy to destroy the regular struc- 
ture of the motion. This point of view is confirmed by the Poincare section 
in oscillator variables which is displayed in fig. 9: For small energies of the 
oscillator the dynamics is regular, whereas in the outer region of the oscillator 
phase plane where its energy is high the trajectories are chaotic. 

Resonance situation 

Another important feature displayed by the full Bloch - oscillator system is 
the possibility of resonance between the excitonic and vibronic subsystems. In 
the DST approximation this effect is absent due to the quasistationarity of 
the vibronic subsystem. In fig. 10 Poincare sections corresponding to such a 
resonance are displayed (r = 1.0, i. e. the transfer frequency equals the vibronic 
frequency 2V = u). It should be noted that the coupling strength p = 0.5 < 1 
is well below the bifurcation value and so there is no hyperbolic point in the 
bonding region. The antibonding fixed point, however, is unstable due to the 
resonance and therefore the dynamics becomes irregular as the energy enters 
the antibonding region. 

4.3 Asymmetric dimer 

The two dimensional character of the bifurcation in the case of an asymmetric 
dimer and the dynamic behaviour of the system associated with it is illustrated 
in the fig. 11. The coupling strength (p = 2) is chosen at a value where the 
symmetric dimer would be above the bifurcaion. The asymmetry parameter 
(g = 0.5) is sufficiently large that the (p,q) pair is outside the bifurcation 
region. Correspondingly the dynamics of the system is regular for intermediate 
energies as in fig. 11a. The effect of the asymmetry is, however, supressed by a 
high oscillator energy as it can be seen in fig. lib where the dynamics for the 
same (p, q) parameter set but at much higher total energy is displayed. There 
is no important difference compared to the symmetric dimer above bifurcation 
from fig. 8c, d. This shows that in the high energy region the dynamic part Q 
is much larger than the static difference a contributing to the shift in the site 
energies (14) and the asymmetry of the coupling strengths q in (18). Hence 
the system is effectively symmetric and behaves irregular because of a coupling 
strength p > 1. 

5 Conclusions 

1. The simplest model of electronic excitation transfer coupled to vibronic 
degrees of freedom, namely the excitation transfer between two sites cou- 
pled to a polarization oscillator possesses a phase space with mixed prop- 
erties: The system displays chaotic or regular behaviour depending on 
the initial conditions on the Bloch sphere and the coupling strength to 
the polarization oscillator. Basic to the properties of the system is the 
bifurcation converting the former ground state in the bonding region of 
the Bloch sphere for a coupling strength above a critical value into an 
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unstable hyperbolic point from which chaotic dynamic starts to develop 
and spreads over the Bloch sphere with increasing total energy of the 
system. The region of antibonding states on the Bloch sphere is the re- 
gion of highest energy in the electronic subsystem where compared to the 
bonding region for a fixed total energy the energy of the oscillator subsys- 
tem is lowest and the regular structures persist destruction until the total 
energy is well above that of the antibonding states. Besides high energy 
and overcritical polarization there is a second mechanism for irregular 
dynamics: resonance between the electronic and vibronic subsystems. 

2. The character of the bifurcation is not changed when asymmetry effects 
in the coupling strength and electronic site energies are introduced. The 
bifurcation diagram becomes, however, two dimensional containing an 
additional asymmetry parameter. The asymmetry parameter can vanish 
when the difference of the coupling strengths is exactly compensated by 
the energy mismatch of the site energies. There is also a dynamic supres- 
sion of asymmetry possible when the changes in the site energies produced 
by the irregular motion of the oscillator in the chaotic region are much 
larger than the differences in the coupling strength and / or site energies. 

3. In the DST approximation the polarization oscillator is assumed to be 
always in its ground state. The type of the bifurcation of the fixed points 
is the same for the full system and the DST case. In the latter case, 
however, the bifurcation is realized on the integrable flow of two indepen- 
dent variables on the Bloch sphere. The hyperbolic point reduces to the 
crossing of the separatrices in a homoclinic structure on the Bloch sphere. 
The perturbation of this homoclinic structure by a non - integrable per- 
turbation is a well known route to chaos. Considering the dynamics of 
the polarization oscillator as such a perturbation the DST solutions can 
be used as a reference to understand the onset of chaotic behaviour in the 
hyperbolic point / separatrix region. 
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